Stability of gold nanowires at large Au-Au separations 
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Abstract 



The unusual structural stability of gold nanowires at large separations of 
gold atoms is explained from first-principles quantum mechanical calcula- 
tions. We show that undetected light atoms, in particular hydrogen, sta- 
bilize the experimentally observed structures, which would be unstable in 
pure gold wires. The enhanced cohesion is due to the partial charge transfer 
from gold to the light atoms. This finding should resolve a long-standing 
controversy between theoretical predictions and experimental observations. 
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Gold is known as the noblest metal of all due to its inert behavior as a bulk material 
and low reactivity of its perfect surfaces. This property made gold probably the most 
valuable and desirable metal in modern human history. At the same time, it still chal- 
lenges scientists with some puzzling features. Among them there are striking structural 
and electrical properties of gold nanocontacts. Such nanocontacts are obtained by the 
mechanically controllable break-junction (MCB) technique or by the scanning tunneling 
microscopy (STM), when the tip is first driven into the gold sample, and then retracted 
forming a nanowire On stretching the nanowire undergoes a number of structural 

rearrangements finally producing a monoatomic gold bridge suspended between the tip 
and sample. In most experiments such nanowires are found to be amazingly stable, reach- 
ing the interatomic distances of about 3.6 - 4.0 A f||-§J i.e. significantly larger than the 
nearest-neighbor distance in bulk gold (2.88 A). The conductance of the nanocontacts 
is quantized in steps of 2e 2 /h approaching the theoretical one-dimensional limit of 2e 2 /h 
when the suspended wire becomes monoatomic A different approach to produce and 
study gold nanowires is based on the transmission electron microscopy (TEM). In this 
technique gold nanowires are fabricated by perforating holes in a gold thin film using a fo- 
cused electron beam with current density ~ 100 A/ cm 2 . The structural reorganization of 
the nanowire is further studied by the electron beam with current density either reduced 
to ~ 30 A /cm 2 or kept the same as in the preparation process ||. The results of 
the former experiments HJ3J] generally agree with the results of the STM and MCB based 
experiments, finding long gold interatomic distances of 3.4 - 3.6 A |4j]. On the contrary, 
the latter experiment || does not confirm the existence of long Au-Au bonds, reporting 
the bulk-like spacing of ~ 2.9 A. 

The issue of the novel gold nanowires has already been addressed in a large body 
of theoretical and experimental works [§^0 which describe many important features of 
gold nanocontacts, including their formation, structural rearrangements and conducting 
properties. However, a convincing explanation of the puzzling structural stability of gold 
monoatomic wires at large Au - Au separations, seen in most experiments, is still lacking. 
In the present report we show from first-principles quantum mechanical calculations that 
the unusual structural stability of gold nanowires at large separations of gold atoms can be 
explained by the presence of undetected light atoms able to stabilize the experimentally 
observed structures, which would be unstable in pure gold wires. In particular, we find 
hydrogen to be a primary suspect. 

To understand the origin of the phenomena we have performed first-principles calcula- 
tions in the framework of the density functional theory (DFT) [[U|. Two sets of data have 
been calculated at T = K, first, for the gold nanowires supported by gold tips consisted 
of 28 atoms, with monoatomic part consisted of 2 to 6 atoms, and second, for infinite 
monoatomic chains with the similar number of gold atoms treated as inequivalent in the 
periodically reproduced unit cells. In agreement with previous findings [pUID], the results 
obtained for both sets (with and without supporting tips) are very similar, indicating 
that the stability of a nanocontact is essentially determined by its monoatomic part. No 
magnetization is found for any studied configuration. The binding energy curve for an 
infinite gold monoatomic wire is given in Fig. 1(a) together with the equilibrium atomic 
configurations and the ranges of their stability. Zigzag is an energetically favorable atomic 
configuration. Upon stretching the wire becomes linear at the Au-Au distance of about 
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2.7 A remaining stable only in a narrow interval up to about 2.9 A, when the Peierls dis- 
tortion triggers the dimerization of the wire accompanied by a metal-insulator transition 
H (Fig. 1(a)). Decreasing the Au-Au distance below the range shown in Fig. 1(a) leads 
to a series of deepening local energy minima and corresponds to the chain folding into the 
bulk-like structure. 

These results suggest the wire to break (and obviously stop being conductive) already 
at interatomic distances of 2.9 A, in agreement with previous calculations [|T| p~0j| and the 
experiment by Takai et. al f|. They, however, by no means explain the stability of gold 
monoatomic wires with the spacing of 3.4 - 4.0 A observed in most of other experiments 
0-^]. Therefore the existing first-principles theory does not provide a complete description 
of the phenomenon. It is especially remarkable regarding the fact that very different gold 
systems (bulk, surface, dimers, atomic clusters), and in particular their geometry, are 
known to be accurately described in the framework of the DFT ||14|| . It seems to be a 
reasonable question to ask if there could be anything missing in the previous theoretical 
considerations of gold nanowires. For instance, it is tempting to check whether small 
amounts of light atoms, which are hard to detect in standard experiments, may influence 
the stability of gold nanocontacts. 

We start with an observation that very similar results for the Au-Au spacing (3.4 - 4.0 
A) have been obtained in different experiments, with and without (ultra) high vacuum 
conditions, in situ sample preparation, and controlled low temperature P-44],|T5|. Espe- 



cially in the case of the ultra-high vacuum the probability to encounter most of the light 
atoms, such as, for instance, carbon and oxygen, is expected to be low. In particular 
it has been shown by Ugarte fl6j and Rodrigues and Ugarte that in the TEM-based 
experiments intense electron irradiation is helpful against carbon contamination. Never- 
theless, such experiments |4j] still show large Au-Au distances in gold nanowires. There is, 
however, at least one element, hydrogen, which is present even at high vacuum conditions 
fl7yi8|l . It is also known that gold becomes unrecognizably active when deposited as small 
particles on a substrate Hl9|-[2l|1. This enhancement of gold reactivity is caused by a high 
concentration of steps and surface defects in small particles, as well as by the substrate- 



induced strain f22]j . Note that such low-dimensional gold structures are well presented on 
the surfaces of the STM tips. Gold surfaces interact with atomic hydrogen, which is able 
to form a stable layer on the surface |[23|| . Dissolving into the bulk, hydrogen tends to 



segregate at the surface and subsurface layers |M] . Gold has also been reported to play a 
role of a hydrogen trap in gold alloys that may be considered as an indication of a strong 
interaction between gold and hydrogen atoms . 

We have extensively studied this interaction from first-principles in the systems of 
different scales: from nanostructures to the fee bulk gold pEfl . Our results show that 
hydrogen molecules do dissociate when embedded into the gold matrix. Moreover, being 
dissolved in small quantities in the fee bulk gold, hydrogen binds strongly to Au atoms 
with characteristic Au-H distances of about 1.75 A. Moderate stress and deformations of 
the gold host lattice do not break these bonds making hydrogen follow Au atoms. 

The most profound influence of the hydrogen presence has been observed in gold 
nanowires. The binding energy curve for a one- dimensional gold-hydrogen wire is pre- 
sented in Fig. 1(b). There is a sequence of structural transformations similar to that of a 
pure gold wire (see Fig. 1(a)). A zigzag arrangement of atoms is most stable for Au-Au 
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distances smaller than 2.7 A, followed upon stretching by the linear configuration of the 
gold atoms with zigzag-like positions of the hydrogen atoms. This configuration is stable 
up to about 3.4 A of Au-Au distance. On further stretching the chain becomes linear 
with hydrogen placed in between the gold atoms reaching the Au-Au spacing of about 
3.8 A and Au - H spacing of 1.85 A. Further stretching results in dimerization eventually 
followed by the complete break of the wire. Thus monoatomic gold-hydrogen wire shows 
an enhanced stability at large Au-Au distances that well matches the experimentally ob- 
served spacing of 3.6 - 4.0 A The cohesion is twice as strong as the one of a pure 
gold chain (see Fig. 1). As experiments on gold nanowires are often performed at room 
temperature, we have also studied the stability of monoatomic Au-H wires at finite tem- 



peratures by first-principles molecular dynamics (MD) simulations [57J. According to our 
MD simulations the structures of the nanowires obtained at T = K remain intact at 
least up to T = 300 K. 

The nature of bonding in gold nanocontacts containing hydrogen can be understood 
by analyzing the charge density as well as the difference between the actual charge density 
and the atomic densities of gold and hydrogen constituting the system (Fig. 2). Charge 
density presented in Fig. 2(a) reveals pronounced electron density bridges caused by 
hydrogen situated in between the gold atoms. This resembles the hydrogen-type bonding 
where the charge transfer from Au to H, which can be seen in the charge density difference 
plot (Fig. 2(b)), is a key component. Therefore, an enhanced effective binding between 
the gold atoms via Au-to-H charge transfer mechanism results in the striking stability of 
the gold-hydrogen nanowires. 

Further, in experiments 0|| the conductance of gold nanocontacts has been measured 
at a low d.c. bias (about 10 mV). We have also studied the stability and ability of the 
gold-hydrogen chains to conduct at a similar bias, which is simulated by shifting the 
chemical potential at one end of the wire with respect to the other in the way described 
in Ref. p8| . The band structure of such a wire (Fig. 3) clearly shows the existence of 
one conduction channel, in agreement with theory and most experimental findings @-@]. 
Additional forces appearing on atoms due to the applied bias are of about 0.001 eV/A 
that is by no means sufficient for breaking the bonds. 

In general one cannot rule out the possibility of the presence of other light elements. 
First, small amount of light atoms can be already present in the sample used to fabricate 
a nanowire. Second, not all experiments on gold nanowires employ high vacuum condi- 
tions. We have considered carbon and oxygen as possible contaminants (see Ref. [|1^] for 
calculational details) and found that both of them stabilize monoatomic wires up to the 
Au-Au distances of ~ 4.5 A. The enhanced cohesion is due to the bonding mechanism 
similar to the one operating in Au-H wires. The equilibrium Au-Au distances in Au-0 
and Au-C wires are, however, somewhat larger (~ 3.6 A in both systems). Thus we notice 
that due to the larger equilibrium Au-Au separation in such systems, they should be in 
a " compressed" state in the region of experimentally observed Au-Au distances, contrary 
to the Au - H system. Therefore, if in a particular experiment oxygen or carbon atoms 
are responsible for the increased stability of gold nanowires, one would at least expect a 
different dependence of the force on the tips displacement from what is actually seen in 
the STM-based experiments |TU| . 



Further, we comment on the experimental results by Takai et al. M reporting the 
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bulk-like spacing in gold monoatomic wires and, therefore, confirming the theoretical re- 
sults for pure gold wires. We suggest that this experiment was most probably free of 
contamination due to a permanent use of intense electron irradiation (current density ~ 
100 A/ cm?) during the observations. As a result, gold nanowires were stable at a bulk-like 
atomic separation of 2.9 A within a very short time interval. It is in a sharp contrast with 
other experiments 0,|J, where an electron beam with current density of 100 A /cm 2 was 
only used to perforate neighboring holes in gold films while the observations were per- 
formed at essentially lower current density (~ 30 Aj cm 2 ). As we already mentioned, these 
experiments resulted in long-living monoatomic wires with large interatomic distances. 

To summarize, we have shown that the striking stability of the suspended gold 
nanowires may be attributed to the hydrogen atoms, which mediate gold-gold interac- 
tions, dramatically increasing binding. Hydrogen, being always present but rarely con- 
trolled in experiments, alters the properties of gold nanowires, the crucial components of 
nanometer-scale devices developed for microelectronics and biophysics applications 



In addition we note that other light atoms, like carbon and oxygen, are also able to sta- 
bilize gold nanowires. Our calculations show that Au-0 and Au-C monoatomic wires are 
stable up to the Au-Au distances of 4.5 A. We suggest that further experiments, incor- 
porating hydrogen and other light atoms detection, may lead to important implications 
for the modern atomic-size technologies. 
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FIGURES 



FIG. 1. Binding energy curves for infinite monoatomic wires of pure gold (a) and gold with 
hydrogen (b). The shown energies are given per gold atom (a) and gold- hydrogen pair (b) as 
a function of the average Au-Au distances. These are calculated as (di +...+ d n )/n, with d n 
being the actual distances between gold atoms in linear chains or the z-projected distances for 
zigzag-like chains H, where n is the number of gold atoms per unit cell. For infinite monoatomic 
wires the total amount of hydrogen atoms (case (b)) is equal to the total amount of gold atoms 
resulting in 2 to 6 hydrogen atoms per simulated nanowire (see Fig. 2). The energy curves for 
the chains forced to be strictly linear (with all d n equal) are shown in blue. The equilibrium 
atomic configurations are given in the corresponding ranges of their stability. The underlined 
intervals indicate the range of stability of the chains with the linear arrangements of gold atoms. 

FIG. 2. Charge density distribution (a) and the difference between the actual density and 
atomic charge densities of gold and hydrogen constituting the system (b) for the monoatomic 
part of the gold-hydrogen nanocontact with the Au-Au distance of 3.6 A. Charge density in (a) 
spans the range 0-1 electron/ A 3 with the spacing of 0.1 electron/ A 3 . Charge density difference 
in (b) spans the range -0.1 - 0.2 electron/A 3 with the spacing of 0.03 electron/A 3 . In both (a) 
and (b), the color background changing from blue via green and red to yellow indicates the 
increase of the density. Accordingly, in the charge density difference plot (b) the blue part of 
the spectrum indicates the lack and the red and yellow ones the excess of electrons. Positions 
of Au and H nuclei are marked with corresponding symbols. 

FIG. 3. The band structure of the gold-hydrogen nanowire with the Au-Au distance of 
3.6 A (see Fig. 1) calculated on a sufficiently dense mesh of k-points (30 k-points in both, TZ 
and TX directions). The energy is given with respect to the Fermi level. The perfectly flat 
bands in TX direction indicate the absence of interaction between single nanowires simulated 
in three-dimensional periodic supercells. 
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